Exact spin-cluster ground states in a mixed diamond chain 



o 
o 



Ken'ichi Takano 1 , Hidenori Suzuki 1 , and Kazuo Hida 2 
1 Toyota Technological Institute, Tenpaku-ku, Nagoya 468-8511, Japan 
2 Division of Material Science, Graduate School of Science and Engineering, 
Saitama University, Saitama, Saitama 338-8570, Japan 
(Dated: 2009.3.31) 

The mixed diamond chain is a frustrated Heisenberg chain composed of successive diamond- 
shaped units with two kinds of spins of magnitudes S and S/2 (S: integer). Ratio A of two exchange 
parameters controls the strength of frustration. With varying A, the Haldane state and several spin 
cluster states appear as the ground state. A spin cluster state is a tensor product of exact local 
eigenstates of cluster spins. We prove that a spin cluster state is the ground state in a finite interval 
of A. For S — 1, we numerically determine the total phase diagram consisting of five phases. 

PACS numbers: 75.10.Jm, 75.10.Pq, 64.70.Tg, 75.30.Kz 
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Effects of frustration under strong quantum fluctuation 
in low-dimensional quantum spin systems have been of 
great interest in condensed matter physics. Among them, 
the diamond chain consisting of successive diamond- 
shaped units has been attracting remarked attention. 

It has been rigorously shown that a ground state of the 
uniform spin diamond chain (UDC) is a singlet spin clus- 
ter state with spontaneous breakdown of translational 
symmetry [l|. The spin cluster state is an essentially 
quantum state described as a tensor product of local sin- 
glet states each consisting of several spins. For a usual 
gapped spin liquid, such a state is an approximation or 
is exact only at a single parameter value. However, in 
the UDC, the spins in each cluster forms an exact local 
singlet, and the exact solution is realized over a finite 
parameter range. 

The distorted version of this model has been also in- 
vestigated theoretically @, H, 0L [f|. Since its discovery 
in a natural mineral azurite [g, [7| , this system is attract- 
ing a renewed interest. Also materials in a unfrustrated 
parameter region are found [1, Q . 

In this paper, we report rigorous ground states includ- 
ing spin clusters which are not singlet. We found the 
cluster states for the mixed diamond chain (MDC) which 
has the same lattice structure as the UDC as shown in 
Fig. [Ha). Unlike the UDC, the MDC consists of two 
kinds of spins with magnitudes S and (S: a pos- 
itive integer). The ratio of two exchange constants in 
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FIG. 1: (a) The MDC, where spin Si has a magnitude S (a 



positive integer) and spin 



_(<*) 



(q = 1,2) has t = |S. The 



solid (dashed) line represents exchange parameter J (K) . (b) 
The equivalent lattice to the lattice in (a) for each set of {Ti}, 
where Ti is the magnitude of Ti = t 



this model is a parameter representing the strength of 
the frustration. For a weak frustration case, the ground 
state of the MDC is the Haldane state , while that of 
the UDC is ferrimagnetic. An interest is how the ground 
state changes from the Haldane state when frustration 
increases. We show that the MDC undergoes a series of 
characteristic quantum phase transitions with increasing 
frustration. A typical ground state is of a tensor prod- 
uct form of finite-length triplet clusters and dimers in 
an alternating array; thus the translational symmetry is 
spontaneously broken. Rigorous results on the phase dia- 
gram are presented for general values of S. In particular, 
a tensor-product state is proven to be the ground state 
for a finite parameter range. For S = 1, the whole phase 
diagram is numerically determined. 

The cluster ground states of the MDC have a macro- 
scopic degeneracy, because each non-singlet cluster de- 
generates with respect to the spin direction. Such macro- 
scopic degeneracy can give birth to a rich variety of quan- 
tum phenomena in the presence of perturbation. Even 
within our preliminary studies, quantum ferrimagnetism, 
unusual Haldane phases with spontaneous breakdown of 
translational symmetry, and infinite series of quantum 
phase transitions are found in the presence of lattice dis- 
tortion. 

The MDC consisting of N diamond units is described 
by the Hamiltonian: 



N 



H 



N 



5> 



1=1 



H t = J(S, + S m ) • (if ) + O + Kt 
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In the ith unit cell, Si is a spin whose magnitude S is an 
integer, while t+ (a = 1, 2) is a spin of magnitude r = 
7jS. Below we sometimes use the notations, Sj = S, + 

S i+ i and = if ) + rf ] . We assume that J > 0. Then 
the ratio A = K/J represents the strength of frustration 
if K > 0. 

We first examine the classical limit of S — > oo with 
finite JS 2 . (i) For A < 2, we write the classical version 
of Eq. © as Uf = \ J[(2T, + S 4 ) 2 - 2(2 - A)T 2 - S 2 - 
AS 2 ]. This is minimized if |T,| = 2r = S, ISA = 2S 
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and \Ti + §S,-| = 0. Then all the Si's (t^'s) in the 
chain are aligned parallel (antiparallel) to a fixed axis, 
and the ground state is antiferromagnctic. This ground 
state is elastic, since any local modification of the spin 
configuration increases the energy, (ii) For A > 2, we use 
the expression Uf = (J/4A)[2(AT< + Si) 2 - 2Sf - A 2 5 2 ]. 
This is minimized if |S,| = 25 and |Ti+Si/A| = 0. Hence 
Si and Si+i are parallel, and t[ and t$ form a triangle 

with Si/A. and r| 2 ' may be rotated about the axis 
of Si and Si + i without raising the energy. Then all the 
Si's in the chain are aligned parallel to a fixed axis, and 
the arbitrariness of the local rotation of t\ and r l - 2 ' ) is 
not obstructed. Thus the ground state is ferrimagnetic 
with magnetization (1 — 2/X)SN. 

Returning to the quantum case with general 5, the 
MDC has a series of conservation laws 

[TlH N } = (i = l,2,...,JV); (3) 

i. e. T{ (= 0, 1, • • • , 5) of T 2 = T^Ti + 1) for each i 
is a good quantum number. With a fixed set {Ti}, the 
original problem of 3N spins reduces to a problem of a 
linear chain with 2N spins, where the (2i — l)th site is 
occupied by the spin Si and the 2zth by Ti. The energy 
of A'-bonds is determined solely by {Ti}. In particular, 
if Ti = 0, then there is no interaction between the left 
and right sides of Ti and the whole lattice is decoupled. 
We denote the lowest energy of Hm for a fixed set {Ti} 
as -E/v({Ti}, A), and sometimes drop {Ti} and/or A in it. 
Hereafter we use an energy unit of J = 1 . 

There is an equivalent lattice to the MDC for each 
fixed set of {Ti}, because T 2 s are conserved. The equiv- 
alent lattice for set {T} is a nearest-neighbor antiferro- 
magnetic linear spin chain on 2N + 1 sites. The spin 
magnitudes are 5 on all odd sites and Ti on 2?th site for 
all i, as seen in Fig. W[b). The Hamiltonian is 

N 

H N = J2(Si + Si+i) • Ti. (4) 

i=l 

H N is related to H N as H N = H N + |A EjliPH^R 1 ) - 
+ 2)]. The ground state energy of H N is denoted 
by E N ({Ti}), or simply E N . 

We begin with analyzing eigenstates of an isolated di- 
amond unit described by Hi- This is realized in the total 
system, if Ti_i = T i+ i = 0. Using the expression, 

7U = i{(T, + S,) 2 + (A - 1)T 2 - S 2 } - ^5(5 + 2), 

(5) 

the lowest energy of Hi for a given Ti is 

^(Ti) = ^(Ti + 1)[A - A(Ti)] - -^5(5 + 2), (6) 

where A(L) = 2(25 + 1)/(L + 1) for nonnegative integer 
L. Then the total spin of the diamond unit is 25 — Tj. 
For the whole spin chain, we have the following lemma: 



Lemma 1. If L > 1 and A > A(L), then T t ^ L for 
any i in the ground state of the MDC. 

Proof. The total Hamiltonian is divided as H = 
H m + H') where N of Hn is omitted and H' is the sum 
of Wi's with i ^ m. We take a state |* ) = |0 m ) ® |*') 
where |0 m ) is the singlet wave function of T m and 
is the lowest-energy state of the other spins. Since 
(0 m \H m \0m) = Ei(0) does not involve S m and S m +i, we 
have (* |-ff|*o) = Ei(0)+E' witbE' = (■$f'\H'\W) being 
the ground-state energy of H' . Let |W) be any state with 
T m = L> 1. Then = (*|W m |*) + (*\H'\iJ/). 

Clearly (^|H m |*) > E\{L) and {V\H'\V) > E'; also 
E X (L) > E^Q) for A > A(L). Therefore (*|-ff|*) > 
(* |-H"l*o) and |*) is not the ground state of H. □ 

Since A(L) decreases with L from A(l) = 25 + 1 to 
A(5) = 2(25+ l)/(5+ 1), we obtain the following result: 

Theorem 1. If A > 25+ 1, then T, = for all i in the 
ground state. 

In this ground state, all pairs of and form 
singlet dimers, so that all Si's are decoupled from each 
other and behave as free spins. Hence we call this state 
the dimer-monomer (DM) state. The DM state for 5=1 
is shown in Fig. HJa). The picture is also the same for 
the DM state of a system with 5 > 2. The ground state 
energy is A^ £i (0) for the lattice consisting of A^ diamond 
units, where 2Af spins with magnitude ^5 and A^ + 1 
spins with magnitude 5 are included. Due to the free 
spins, there is a (25 + l) Ar -fold degeneracy in the DM 
state. 

For A < 25 + 1, we have the following lemma: 

Lemma 2. For A in §(25 + 1) < A < 25 + 1, T t at 
any site is or 1 in the ground state. Also at least one 
of TiS is not 0. 

Proof. For r = |, the ground state is composed of 
TiS with their magnitude or 1. It is true even for r > 1 
if A(2) < A < A(l) because of Lemma 1; the region 
is rewritten as |(25 +1) < A < 25 + 1. Further for 
A < 25 + 1, the DM state, where Ti = for all i, is 
not a ground state, because we can lower the energy by 
introducing an isolated Ti = 1 in the DM state as known 
from Eq. ©. □ 

Hence Aci = A(l) = 25+1 is the phase boundary of 
the DM phase. 

For | (25 + 1) < A < 25 + 1, the ground state consists 
of an appropriate array of Ti = and 1. For further 
arguments, we define the cluster-n as the local ground 
state of a successive n diamond units in the case that T;S 
are all unity in the n diamond units and both TiS just 
outsides of it are zero; the cluster-n is a state of n+ 1 SiS, 
n rj^s and n T^s. The magnetization of the cluster- 
n is (n + 1)5 — n due to the Lieb-Mattis theorem [TT[ | . 
In particular, for 5=1, the cluster-n is the Haldane 
state of a finite Haldane chain with 2n + 1 sites, and is 
a triplet state. Denoting the energy of the cluster-n by 
E n , we define the energy per diamond unit for a cluster-n 
measured from the DM state as 

e„(A) = — J— [E n -nE 1 (0)]. (7) 
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FIG. 2: The ground states for the MDC in the case of S = 1 
and t = |: (a) The dimer-monomer (DM) state. An un- 
shaded oval represents a singlet dimer (T = 0). There are 
free spins (monomers) on the sites not enclosed by ovals, (b) 
The dimer-cluster-1 (DC1) state. A shaded oval represents 
a triplet pair (Ti — 1). A loop including a diamond unit is 
a cluster- 1 whose spin magnitude is unity, (c) The dimer- 
cluster-2 (DC2) state. A loop including two diamond unit is 
a cluster-2 whose spin magnitude is unity, (d) The dimer- 
cluster-3 (DC3) state. A loop including a diamond unit is a 
cluster-3 whose spin magnitude is unity, (e) The Haldane (H) 
state, (f) The ground-state phase diagram. 



We have divided it by n + 1 and not by n, since a cluster- 
n is accompanied with the neighboring Ti = 0s, and it 
is convenient to include the effect of one of them. If 
e„(A) is the minimum at a single value of n, then the 
ground state in the thermodynamic limit (N — > oo) is 
realized by an alternating array of cluster-ns with Tj = 1 
and dimers with Ti = 0. We call the state the dimer- 
cluster-n (DCn) state. The DCn states with n — 1, 2, 
and 3 for 5 = 1 are illustrated in Fig. [2|b), (c), and (d), 
respectively. The DCn state has spatial periodicity n+1, 
and the translational symmetry is spontaneously broken. 
The degeneracy is (n + l)3 N '^ n+1 ' for large N. 

To investigate the ground state of the MDC in the 
region of T s ; = or 1 for all i, we need the following 
lemma for small equivalent lattices: 

Lemma 3. E 2 > 2E U if T. L = 1 for all i. 

Proof. We define three spin Hamiltonians H\l = 
(Si + S a ) • Ti and H 1R = (S 2 + S 3 ) • T 2 . The five 
spin Hamiltonian for Si, Ti, S2, T 2 , and S3 is writ- 
ten as H 2 = H\l® It 2 ,s 3 + ir^Si ®Hir, where It,s is an 
identity operator in the subspace of spins T and S. The 
ground states of H\l, H\ R , and H2 are denoted by \4>3l), 
|^>3.r), and IV's)) respectively. From variational consider- 



ation, the following inequality holds: E2 = (V'sl-^lV's) > 
(i>3L\H 1L \iP 3L ) + (*p 3R \H 1R \ip 3R ) = 2Ei. 

We now deduce a contradiction by assuming E2 = 2E\ . 
This is valid only if ^5) is also the ground state of both 
Hit, ® It 2 ,s 3 arid It u Si ®H\ R . The Lieb-Mattis theorem 
[111 ] applied to H\l says that the magnitude of the total 
spin Si + Ti + S 2 must be 25 - 1 in \ip 3L )- Therefore 
|^ 3i ) contains only the states with Sf+T{+S$ < 25-1. 
Applying the Lieb-Mattis theorem to H%, however, the 
ground state ^5) has total spin 35 — 2 and contains all 
the 5 z -diagonal states \Sf,Tf, 5f , T|, 5f ) with total spin 
35 — 2. This implies that ^5) contains the state with 
Sf + Tf + SI = 25 such as |5, 0, 5, -1, 5 - 1) with finite 
amplitude. This is a contradiction and therefore we have 
E 2 > 2Ei. □ 

Based on Lemma 3, we have the following proposition 
on the ground state in a region just below Aci = 25+1: 

Theorem 2. There exists a positive number S such 
that the DC1 state is the ground state for Aci — S < A < 
Act in the thermodynamic limit. 

Proof. Due to Lemma 2, we only have spin magni- 
tudes or 1 for T,s for §(25 + 1) < A < 25 + 1. We 
consider a segment of n diamond units between two 0s. 
The equivalent lattice is described by Eq. ((4]) with N = n 
and Ti — 1 for all i. By using the ground state energy E n , 
we have E n — nE% (0) = E n + n\. We adopt 5 = E 2 — 2Ei 
(> 0) if it is smaller than |(25+1), and show that d sat- 
isfies the condition of Theorem 2. 

If n is even, the Hamiltonian for E n is divided into 
^ sub-Hamiltonians, each equivalent to that for E 2 . A 
variational argument gives E n > ri_B 2 /2 = nEi+nd/2 for 
even n, and E n > (n - l)E 2 /2 + Ei = nE x + (n - 1)5/2 
for odd n. Since e% = (E\ + A)/2 = (A — Aci)/2, a lower 
bound on e n — e% is given as 2 (ri+i) — wn i cn 
is positive for A > Aci — 6 and n > 2. 

If 5 = E 2 - 2^i > |(25 + 1), we use 5' = |(25 + 1) 
instead of 6. Then Aci — 6' is the lower bound of Lemma 
2 where spin magnitudes of T^s are or 1. The argument 
of the last paragraph stand still for 6' replacing 8 since 
< 5' < 5. □ 

For A < 1, we show that the ground state of the MDC 
is equivalent to that of the uniform linear spin chain with 
spin magnitude 5. 

Theorem 3. If A < 1, then T,- = 5 for all i in the 
ground state. 

Proof. We divide the total lattice into two sublattices 
A and B, where SiS are on the A sublattice and t\ s 
on the B sublattice. For A < 1, we have K < J for 
interaction J (> 0) between spins on the different sub- 
lattices and interaction K between spins on the same B 
sublattice. This suits to the condition for the Lieb-Mattis 
theorem [ll[ . Hence the total spin of the ground state is 
given as \(N + 1)5 — NS\ = 5 for the lattice consisting 
of diamond units. 

We also apply the Lieb-Mattis theorem to the ground 
state of the equivalent lattice with {Ti}. Then the total 
spin of the ground state is given as | (N + 1)5 — J2i=i Ti | • 
This must be the same value 5 as the total spin of the 
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FIG. 3: A-dependence of the (a) Curie constant and (b) resid- 
ual entropy for the MDC with S — 1. 



ground state of the MDC. This is possible only if T, = S 
for all i. □ 

Thus, at least for A < 1, the MDC is equivalent to a 
uniform linear chain with integer spin magnitude S. The 
gro und state of an integer spin chain is the Haldane state 
(lfj | . The picture of the Haldane state for the MDC with 
S = 1 is shown in Fig. [2je). 

For S — 1, the DCn ground state is resolved into the 
ground states of spin clusters equivalent to finite-length 
spin-1 chains. Therefore we can precisely determine the 
phase boundaries by the exact numerical diagonalization 
for the finite chains. Equation (O becomes 

where E n is the ground state energy of the spin-1 chain 
with length 2n + 1. Typical values are E$ = 0, E\ = —3, 
E 2 = -5.8302125227708, and E 3 = -8.6345319827062. 
If Ti = 1 for all i, the ground state energy per unit cell is 
given by e«,(A) = 2e oa + A, where e m ~ -1.401484038971 
[12j | is the ground state energy of an infinite spin-1 
chain per unit cell. The phase transition between the 
DC(n — 1) and DCn phases takes place at A = Acn = 
(n + l)E n -i —nE n , which is the solution of e„_i(Acn) = 
en(Acn)- However the DCn phases with n > 4 do not 
appear [l3T | and a phase transition takes place directly 
from the DC3 phase to the Haldane (DCoo) phase at 
Acoo = E$ — 8eoo, which is the solution of e3(Acoo) = 



eoo(Acoo)- The critical values of A are estimated as Act = 
25+1 = 3, A C2 = 2.660425045542, A C3 = 2.58274585704, 
and Acoo = 2.5773403291. Thus we have the ground- 
state phase diagram for S = 1 in Fig. Hff). Magneti- 
cally, each phase is characterized by the number of alive 
spins; it equals the number N/(n + l) of triplet clusters in 
the DCn ground state. Consequently, as an experimen- 
tally measurable quantity, the Curie constant given by 
C = |(n+ 1) _1 shows a stepwise A-dependence as shown 
in Fig. 02a). Also, the residual entropy So per unit length 
is ln3/(n + 1) in the DCn phase. At the phase bound- 
ary between the DCn and DC(n -I- 1) phases, the mixing 
of cluster-n and cluster- (n + 1) remarkably enhances the 
entropy, which is estimated by combinatory argument as 
1.333, 0.744, and 0.522 for n = 0, 1, and 2, respectively. 
Thus So na s spike-like structures as shown in Fig. \5j[b). 
Similar structures are also reported in the diamond hier- 
archical Ising model [3| • 

Summarizing, we investigated the MDC and proved 
that the DM, DC1, and Haldane phases exist for an ar- 
bitrary integer spin magnitude S. For S = 1, we nu- 
merically determined all the phase boundaries of the full 
phase diagram, which consists of the DM, DC1, DC2, 
DC3, and Haldane phases. Since each spin cluster for 
magnitude S has 25+1 degrees of freedom, the total 
ground state is massively degenerate. If this degener- 
acy is lifted by an appropriate perturbation, a variety of 
exotic phases are expected. 
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